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Abstract 

The structural and elastic properties as well as phonon-dispersion relations of the Heusler-type 
alloys Fe2VAl and Fe2VGa are computed using density-functional and density-functional pertur- 
bation theory within the generalized-gradient approximation. The calculated equilibrium lattice 
constants agree well with the experimental values. The elastic constants of Fe2VAl and Fe2VGa 
are predicted for the first time. From the elastic constants the shear modulus, Young's modulus, 
Poisson's ratio, sound velocities and Debye temperatures are obtained. By analyzing the ratio 
between the bulk and shear modulii, we conclude that both Fe2VAl and Fe2VGa are brittle in 
nature. The computed phonon-dispersion relation shows that both compounds are dynamically 
stable in the LI2 structure without any imaginary phonon frequencies. The isomer shifts of Fe in 
the two compounds are discussed in terms of the Fe s partial density of states, which reveal larger 
ionicity/less hybridization in Fe2VGa than in Fe2VAl. For the same reason the Cauchy pressure is 
negative in Fe2VAl but positive in Fe2VGa. 



9 



I. INTRODUCTION 



The Heusler-type intermetallic compounds Fe2VAl and Fe2VGa have recently at- 
tracted great attention due to their facinating thermal, electrical, magnetic and transport 
properties,-^ 1 ^ not only from the basic science prespectives but also from the potential 
application as thermoelectric materials.-^ Though Heusler-type intermetallics commonly 
appear as metals,- semiconducting behaviour has been observed in Fe2VAl and Fe2VGa, as 
evidenced by their negative temperature coefficient of resistivity^. This feature has been 
attributed to the appearance of a pseudo gap in the density of states in the vicinity of the 
Fermi level, and these materials have been characterized as semimetals. Nuclear magnetic 
resonance experiments and optical conductivity measurements further confirmed the exis- 
tence of deep pseudo gaps near the Fermi energy (E^) in both compounds.-^ Recently, it 
has been reported by Lue et al.— that the partial substitution of Ga by Ge in Fe2VGa effec- 
tively dopes electrons to the system thereby leading to a dramatic decrease in the electrical 
resistivity A similar study for Fe2VAl was carried out by Nishino et al.— with the same 
conclusion. The variation of the Seebeck-coefficient with a sign change from possitive to 
negative can be understood by means of a rigid band like shifting of the Fermi level across 
the pseudogap. 

Several band structure calculations have been reported providing various insights into 
Fe2VAl and Fe2VGa, all confirming that Fe2VAl is a non-magnetic semimetal with a pro- 
nounced pseudo gap at the Fermi leve^^^i^yLLL^. The calculated pseudo gap seems 
however too wide to explain the experimental data. Guo et al.— suggested the negative tem- 
perature coefficient of resistance to arise from carrier localisation, while Singh and Mazin^ 
pointed towards the interaction of carriers with fluctuating magnetic moments, and Weht 
and Pickett^ proposed dynamic correlations between holes and electrons to be responsi- 
ble for the unusual resistivity behavior. The importance of magnetic moment formation in 
off-stoichiometric compounds have been confirmed experimentally™^ and recent theoretical 
studies of the Fe2+ x Vi_ x Z (Z=Al,Ga) compounds^ 1 ^ with a supercell approach predicted 
that antisite defects or excess atoms in the Heusler lattice may induce a radical change in 
the electronic and magnetic properties. Similarly, Antonov et al.— have studied the elec- 
tronic structure and X-ray magnetic circular dichroism in Fe2-zVi +a; Al for various x using 
first principles calculations, which confirm the formation of magnetic moment on vanadium 
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impurities. 

Alloys based on Fe2VAl are good candidates for thermoelectric materials.- 1 ^ Although 
they have a high thermoelectric power, because of their high electrical conductivities and 
high Seebeck coefficients, they have poor figures of merit compared with other state-of-the- 
art thermoelectric materials due to their high thermal conductivity, k. Hence reduction 
of k is required for practical applications, 6 which may be achieved by sustitution of Al 
by a heavier element.- Large experimental efforts have been devoted to characterization 
of Fe2VAl based materials,— including x-ray absorption 24 and photoemission 2 ^! to reveal 
their surface and bulk electronic structure and electrical and thermal properties.— Yet, 
the lattice dynamics and mechanical properties have not been explored, which is taken as 
the objective of the present work. Here we present results of computations of the elastic 
constants and phonon- dispersion relations for Fe2VAl and Fe2VGa using density functional 
and density functional pertubation theory within the generalised gradient approximation. 
From the modulus of elasticity, we predict the mechanical behaviour of these compounds. 

II. COMPUTATIONAL DETAILS 

The all-electron linear muffin-tin orbital method- in a full-potential (FP-LMTO) 
implementation^ is used to calculate the total energies, and basic ground state and elastic 
properties of Fe2VAl and Fe2VGa. In this method, the crystal volume is split into two re- 
gions: non-overlapping muffin-tin spheres surrounding each atom and the interstitial region 
between the spheres. We used a double k spdf LMTO basis to describe the valence bands, 
i.e. Hankel-f unctions characterized by decay constants k are smoothly augmented with nu- 
merical radial functions within the atomic spheres. The calculations included the 4s, 4p, 3p, 
and 3d basis functions for iron and vanadium, the 3s, 3p and 3d basis for aluminum, and the 
4s, 4p and 3d bases from gallium. The exchange correlation potential was calculated within 
the generalized gradient approximation (GGA) scheme^. Inside the muffin-tin spheres, the 
charge density and potential were expanded in terms of spherical harmonics up to / ma x=6, 
while in the interstitial region, they were expanded in plane waves, with 12050 waves (en- 
ergy up to 171 Ry) included in the calculation. Total energies were calculated as a function 
of volume for a (28 28 28) k-mesh containing 624 k-points in the irreducible wedge of the 
Brillouin zone. The energy curves were fitted to a second order Birch equation of state 33 to 
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FIG. 1: Calculated phonon dispersion curves and phonon density of states of Fe2VAl 



xom the variation 
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obtain the ground state properties. The elastic constants were obtained 
of the total energy under volume-conserving strains, as outlined in Refs. 

In order to calculate the vibrational properties of the Heusler alloys we used the density 
functional perturbation theory as implemented in the plane-wave pseudopotential method, 
through the Quantum-ESPRESSO packag We have used the Troullier-Martins non- 
conserving pseudopotentials^.. Convergence tests lead to the choice of kinetic energy cutoffs 
of 80 Ry, and a (8 8 8) Monkhorst-Pack 34 grid of k-points for the Brillouin zone integration. 



III. PHONONS 



The calculated phonon dispersion curves of Fe2 VA1 and Fe2VGa are presented in figures 1 
and 2, respectively. We do not find any imaginary phonon frequency in the whole Brillouin 
zone for the two compounds. This supports the dynamical stability of both compounds in 
the Heusler structure, which is not a-priori evident in view of instabilities found in similar 
alloys, e. g. Ni 2 MnAl and Ni 2 MnGa.— ^ It is interesting to note that the optical phonons 
for Fe2VGa are coupled, falling in a broad frequency range between 5.2 and 9.1 THz. This 



is in contrast to Fe2VAl, where three optical phonon branches are separated from the lower 
frequency phonons. Inspection of the atomic masses of Fe, V, Ga and Al allows us to 
understand this difference between the two compounds. The similar atomic masses of Fe, 
V and Ga lead to the coupled phonon dispersions, while the significantly lower atomic mass 
of Al results in the phonon seperation. The three highest optical phonon modes of Fe2VAl 
have their major amplitude contribution on Al atoms. 

The lattice thermal conductivity is inversely proportional to the sound velocity,™ so as- 
suming similar phonon scattering properties (relaxation times) for Fe2VAl and Fe2VGa, the 
former compound thus should have the smaller thermal conductivity, which might be an 
observation of relevance in efforts of minimization of thermal conductivity for good thermo- 
electric figures of merits. Note that in the study of Ge substitution for Al in Fe2VAl in Ref. 
[5 the opposite trend was found, i.e. k is reduced when the heavier Ge substitutes for Al, 
but this was attributed to the change in phonon relaxation time due to the incorporation of 
impurities with a large mass difference. 

IV. GROUND STATE AND ELASTIC PROPERTIES 

The calculated equilibrium lattice constant and the bulk modulus along with the exper- 
imental and other theory work is given Table. 1. The lattice constant obtained from the 
present calculation is underestimated by 0.8% for Fe2VAl and 0.7% for Fe2VGa when com- 
pared with the experimental value. The present results on lattice constant and bulk modulus 
agree quite well with the recent FLAPW calculations. 24 The calculated elastic constants of 
Fe2VAl and Fe2VGa are presented in Table II. This is the first theoretical prediction of the 
elastic constants, which still awaits experimental confirmation. It is noticeable that the 
Cauchy pressure, Cp = Cu — C44 is negative for Fe2VAl (Cp = —45 GPa) but positive for 
Fe2VGa (Cp = +9 GPa). According to Pettifor— this reflects a more covalent character 
of bonding in Fe2VAl. A simple relationship, empirically linking the plastic properties of 
materials with their elastic moduli has been proposed by Pugh.— The shear modulus G 
represents the resistance to plastic deformation, while the bulk modulus B represents the 
resistance to fracture. A high B/G ratio may then be associated with ductility, whereas a 
low value would correspond to a more brittle nature. The critical value separating ductile 



and brittle materials is around 1.75 according to Ref. 
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In the case of Fe2VAl the value 



of B/G is 1.37 and for Fe2VGa the value is 1.73 from our calculated values in Table II, 
and therefore Fe2VAl can be classified as a brittle material, while Fe2VGa resides on the 
borderline between brittle and ductile. Frantsevich^ 3 - has a similar criterion based on the 
Poisson ratio, namely u < 1/3 (u > 1/3) for brittle (ductile) character, hence categorizing 
both compounds as brittle (see Table II). The difference merely reflects the imprecise nature 
of the concept.— 

The sound velocities of Fe2VAl and Fe2VGa may be derived from the calculated elastic 
constants™* see Table III. The properties of the two compounds are quite similar apart 
from the heavier mass of Ga compared to Al, which is the dominating factor leading to 
smaller sound velocities in Fe2VGa compared to Fe 2 VAl. The Debye temperature may be 
estimated from the simple isotropic approximations, ks&D — fiv m ko, where ks and ku 
are the Boltzmann constant and the Debye vector, respectively, and v m the average sound 
velocity. Hence, the smaller sound velocity in Fe2VGa directly leads to a smaller Debye 
temperature. The sound velocities may also be read from the acoustic phonon branches in 
figures 1 and 2. They are listed in the Table as well. One notices quite some anisotropy, the 



sound velocities along [110] being ~ 30 % larger than along the [100] and [111 



directions. 



45. The isomer 



The isomer shift of Fe in Fe2VAl and Fe2VGa have been measured by Ref. 
shift is proportional to the difference in electron contact density, p(0), between the sample 
and reference materials, 

I.S. = a(p(0)-p re/ (0)). 

We calculated the electron contact density of Fe in the Fe2VAl and Fe2VGa compounds 
and compared to that of cn-Fe. Using the calibration constant a = —0.22a 3 ) mm/s (Ref. 



461 ). this leads to the isomer shifts listed in Table IV. The difference in contact density is 
quite small, being larger in Fe2VAl than in Fe2VGa by 0.45 a^ 3 . Thus the isomer shift of 
the latter compound is larger by 0.10 mm/s, in excellent agreement with the experimental 
finding. However, the calculated isomer shifts relative to a-Fe are about 0.04 mm/s smaller 
than measured, which could be a temperature related effect, since the experiments are 
conducted at room temperature, while our theory pertains to K. The origin of the different 
value for the isomer shift in Fe2VAl and Fe2VGa is the difference in electronegativity of Al 
and Ga. Although the electronic structures of these compounds are very similar,— the 
greater electronegativity of Ga compared to Al transcripts into a deeper position of the 
Ga-dominated deepest s-band in Fe2VGa compared to the position of the equivalent Al 
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FIG. 2: Calculated phonon dispersion curves and phonon density of states of Fe2VGa 

s-band in Fe2VAl. This leads to reduced hybridization of Fe s into this band in Fe2VGa 
compared to Fe2"VAl. To illustrate this, we plot in figure [3] the partial density of states of 
Fe s-character for the two compounds. This is the relevant quantity to consider, since the 
dominant contribution to the electron contact density variations in solids stems from the 
valence s-electrons (only s-states - and relativistic pi/ 2 -states extend their amplitude into 
the nuclear region). Thus, in figure [3] the relevant s-band is the the lower part of the density 
of states (between -10.3 eV and -7.5 eV in Fe2VGa, between -9.3 eV and -6 eV in Fe2VAl, 
relative to the Fermi level). The reduced hybridization of Fe s into this band leads to the 
smaller band width in Fe2VGa compared to Fe2VAl, as well as the fewer integrated number 
of Fe s electrons for this band (Table IV), which almost exactly reflects the difference in 
isomer shifts. In contrast, the upper part of the Fe s partial density of states in figure [3l 
which describes the hybridization of Fe s into the Fe and V <i-bands, is much more similar 
in the two cases, and also integrates to almost the same number of Fe s charge. The larger 
hybridization in Fe2VAl compared to Fe2VGa also explains the difference in sign of the 
Poisson ratio, discussed above. 
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FIG. 3: Calculated Fe s partial density of states (in units of electrons per eV and per Fe atom) in 
Fe2VAl (red) and Fe2VGa (blue). The zero of energy is placed at the Fermi level, which falls in 
the pseudo gap within the Fe and V d-bands. 

V. CONCLUSIONS 

Ab initio electronic structure and phonon frequency calculations based on the density fuc- 
ntional theory have been presented for the Heusler alloys Fe2VAl and Fe2VGa. The elastic 
constants have been obtained, while no experiments exist for comparison. The differences 
in phonon dispersions for the two compounds can be explained by the heavier mass of Ga 
compared to Al, while the difference in Fe isomer shift is attributed to the higher electroneg- 
ativity of Ga compared to Al, the ensuing reduced covalency in Fe2VGa compared to Fe2VAl 
also explaining the difference in Cauchy pressure. 
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TABLE I: Calculated lattice constants in A, bulk modulus in GPa, and its pressure derivative B' 
of Fe2VAl and Fe2VGa at the theoretical equilibrium volume. 



Compound 




Lattice Constant 


Bulk Modulus B' 


te2 VA1 


Expt. 


5.76 






Theory, this work 


5.712 


220.8 5.3 




Other theory 


5.712 a 


212 a 


Fe 2 VGa 


Expt. 


5.77 






Theory, this work 


5.727 


228.5 4.9 




Other theory 


5.7_2^ a 


214 a 






a : Reiki; 




TABLE II: 


Calculated elastic constants, 


shear modulus (G), and Young's modulus (E) all expressed 


in GPa, and Poisson's ratio v for Fe2VAl and Fe2VGa at the theoretical equilibrium volume 




Cll C12 




G E v 


Fe 2 VAl 


415.7 125.3 


170.7 


160.5 387.6 0.208 


Fe 2 VGa 


363.1 161.0 


152.3 


131.8 331.7 0.258 



in 



TABLE III: Calculated longitudinal, shear, and average wave velocity (vi, v s , and v m , respectively) 


in km/s, and the Debye temperature, 


#D, in Kelvin from the average elastic wave velocity for Fe2VAl 


and Fe2VGa 


at the theoretical equilibrium volume. Also listed 


are the directional sound velocities 


read from the acoustic phonon branches along 


directions [100], 


[111] and [110]. 


Compound 






v s 




Fe 2 VAl 


isotropic 


8.13 


4.93 


5.45 447 




[100] 


7.49 


4.44 






[111] 


7.72 


4.46 






[110] 


10.23 


6.41 




Fe 2 VGa 


isotropic 


7.09 


4.05 


4.50 369 




[100] 


6.29 


3.52 






[111] 


6.26 


3.13 






[110] 


8.15 


4.52,4.19 





TABLE IV: Calculated and measured isomer shifts of Fe2VAl and Fe2VGa, in mm/s relative to 
a-Fe. Also listed are the partial charges inside the Fe muffin-tin sphere (radius R = 2.305 a.u.). 
The Fe s-charge within the lowest s-like band is listed quoted in paranthesis (see text). 





Fe 2 VAl 




Fe 2 VGa 


I.S. theo 


+0.021 




+0.121 


I.S. expt a 


+0.058(5) 




+0.161(5) 


Fe n s 


0.394(0.140) 


0.370(0.111) 


n p 


0.437 




0.419 


nd 


6.112 




6.121 


n f 


0.015 




0.023 


a : Ref. 
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